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interferometric: polarization 

CONTROL 

RELATED APPLICATIONS 

5 

This application claims the benefit of U.S. Provisional 
Application Ser. No. 60/692,713, filed on Jun. 20, 2005, 
under 35 U.S.C. §1 19(e). 

ORIGIN OF THE INVENTION io 

The invention described herein was made by one or more 
employees of the United States Government and may be 
manufactured and used by or for the Government of the 
United States of America for governmental purposes without 1 5 
the payment of any royalties thereon or therefor. 

FIELD OF THE DISCLOSURE 

This disclosure relates generally to measurement appara- 2 o 
tus, modeling and techniques for polarization characteriza- 
tion or measurement of signals, in particular to a quasioptical 
phase modulator for polarization control and modulation, and 
more particularly, to techniques and apparatus for variable 
delay polarization modulation. 25 

BACKGROUND 

Astronomical polarimetry is an area currently undergoing 
much study and development, at least in part responsive to 30 
high-sensitivity searches for “B-modes” of the cosmic micro- 
wave background radiation polarization. 

The cosmic microwave background is polarized at the level 
of a few microkelvins. There are two types of polarization, 
respectively known as E-modes and B-modes. The relation- 35 
ship between these modes may be analogous to electrostatics, 
in which the electric field (E-field) has a vanishing curl and 
the magnetic field (B -field) has a vanishing divergence. The 
E-modes arise naturally from Thomson scattering in an inho- 
mogeneous plasma. The B-modes, which are thought to have 40 
an amplitude/magnitude of at most a 0. 1 pK, are not produced 
from plasma physics alone. 

Detecting the B-modes is extremely difficult, particularly 
given that the degree of foreground contamination is 
unknown, and the weak gravitational lensing signal mixes the 45 
relatively strong E-mode signal with the B-mode signal. 

B-modes are signals resulting from cosmic inflation and 
are determined by the density of primordial gravitational 
waves. B-modes thus provide signatures for gravitational 
waves associated with the inflationary epoch and are expected 50 
to provide a direct measurement of the energy scale of infla- 
tion. Amplitudes for B-modes are theorized to be on the order 
of 10 “ 7 to 10 -9 of that of the cosmic background radiation, 
and thus measurement of the B-modes requires a robust 
modulation strategy and effective control over systematic 55 
artifacts. 

Emission from magnetically -aligned dust in our Galaxy 
contributes to interference that will have to be understood in 
order to clearly distinguish and extract the B-mode from the 
total signal. However, this polarized emission also provides a 60 
tool for analyzing the role of magnetic fields in star formation. 
The advent of multiple wavelength submillimeter and far- 
infrared photometers, such as SCUBA2 (a new generation 
submillimeter imager for the James Clerk Maxwell Tele- 
scope) and HAWC (a far-infrared camera for the Strata- 65 
spheric Observatory For Infrared Astronomy (SOFIA)), pro- 
vides opportunity to expand such study. Polarization modules 
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have been developed to facilitate leveraging of these new 
photometric tools for such applications by allowing them to 
function as polarimeters. 

Partial polarization results from statistical correlation 
between the electric field components in the plane perpen- 
dicular to the propagation direction. Such correlation is rep- 
resented via complex quantities, and, as a result, in measure- 
ments of polarized light, it is convenient to employ linear 
combinations of these correlations, such as Stokes param- 
eters, e.g., I, Q, U and V. 

The polarization state of radiation through an optical sys- 
tem may be modeled by determining the transformations that 
describe the mapping of the input to the output polarization 
states. In modeling the types of optical elements associated 
with the polarization modules described in this disclosure, 
Stokes I is decoupled from the other Stokes parameters. For 
this class of elements, the polarization P, as described with 
reference to Eq. (1) below, 

p 2 =2 2 +6^+r 2 , (i) 

is constant. Eq. (1) may be interpreted to describe the points 
on the surface of a sphere in three-dimensional space having 
Q, U and V as coordinate axes. This sphere is known as the 
Poincare sphere, and the action of any given ideal polarization 
modulator may be described by a rotation and/or an inversion 
in this space. Such operations correspond to introduction of a 
phase delay between orthogonal polarizations, and that is the 
physical mechanism operative in a polarization modulator. 
The two degrees of freedom of any given transformation are 
the magnitude of the introduced phase delay and a parameter 
describing the basis used to define the phase delay. These two 
parameters directly define the orientation and the magnitude 
of the rotation on the Poincare sphere: the rotation axis is 
defined by the sphere diameter connecting the two polariza- 
tion states between which the phase is introduced, and the 
magnitude of the rotation is equal to that of the introduced 
phase. 

In order to measure the polarized part of a partially -polar- 
ized signal, it is useful to separate the polarized portion of the 
signal from the unpolarized portion. This is especially useful 
when the fractional polarization is small. One way to do this 
is to methodically change, or modulate, the polarized portion 
of the signal (by changing one of the parameters of the polar- 
ization modulator) while leaving the unpolarized portion 
unaffected. Periodic transformations in Poincare space can 
accomplish this encoding of the polarized portion of the sig- 
nal for subsequent demodulation and detection. A convenient 
way of formulating the problem is to envision a detector that 
is sensitive to Stokes Q when projected onto the sky in the 
absence of modulation. The polarization modulator is then 
systematically changing the polarization state to which the 
detector is sensitive. By measuring the output signal, the 
polarization state of the signal or light may be completely 
characterized. 

One conventional way to implement such a polarization 
modulator is by use of a dielectric birefringent plate. A bire- 
fringent plate comprises a piece of birefringent material cut 
so as to delay one linear polarization component relative to 
the other by the desired amount (generally either one-half or 
on quarter of the wavelength of interest). In this case, the 
phase difference is fixed, and the modulation is accomplished 
by physically rotating the birefringent plate (and hence the 
basis of the introduced phase). 

However, a birefringent plate may be built to measure 
either circular or linear polarization, but cannot measure both. 
Additionally, polarization modulators built using this 
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approach are not readily retuned for use at multiple wave- 
lengths. Further, the requirement to be able to rotate the 
birefringent plate engenders need for a complex ensemble of 
shafts, bearings and gears. 

For the reasons stated above, and for other reasons dis- 5 
cussed below, which will become apparent to those skilled in 
the art upon reading and understanding the present disclosure, 
there are needs in the art to provide improved phase modula- 
tors in support of increasingly stringent and exacting perfor- 
mance and measurement standards in settings such as astro- 10 
nomical observation. 

SUMMARY 

The above-mentioned shortcomings, disadvantages and 15 
problems are addressed herein, which will be understood by 
reading and studying the following disclosure. 

In one aspect, the disclosure encompasses a signal condi- 
tioning module configured for insertion into a photometric 2Q 
system to realize a polarimeter. The module may include a 
cascaded series of variable delay polarization modulators. 
Each modulator of the series includes may include an input 
port, and a first arm formed to include a first reflector and first 
rooftop mirror arranged in opposed relationship. The first 25 
reflector can direct an input radiation signal from the input 
port towards the first rooftop mirror. Each modulator may also 
include an output port and a second arm formed to include a 
second reflector and second rooftop mirror arranged in 
opposed relationship. The second reflector can guide a signal 3Q 
from the second rooftop mirror towards the output port to 
provide an output radiation signal. A beamsplitting grid may 
be placed between the first reflector and the first rooftop 
mirror, and may also be placed between the second reflector 
and the second rooftop mirror. A translation apparatus can 35 
provide adjustment relative optical path length vis-a-vis the 
first arm, the second arm and the grid. 

In another aspect, a millimeter wave receiver may include 
an input antenna, and a polarization modulation chain 
coupled to the input antenna. The modulation chain can 40 
include multiple cascaded variable-delay polarization modu- 
lators, each having an input realizing a fixed basis of phase 
introduction, and a translation apparatus providing variable 
magnitude of phase delay between two orthogonal linear 
polarizations of an input signal. The modulation chain may 45 
further include a beamsplitter, a first arm and a second arm 
optically coupled to one another via the beamsplitter. Each of 
the first and second arms may include opposed reflective 
elements disposed on opposite sides of the beamsplitter. The 
opposed reflective elements in each arm may have at least one 50 
mirror coupled to the translation apparatus such that a first 
distance separating the beamsplitter from at least one of the 
opposed reflective elements is variable. 

In a further aspect, the present disclosure contemplates a 
radio telescope, which may comprise an antenna and a 55 
receiver coupled to the antenna. The receiver can provide a 
polarizing interferometer function by including a polarization 
modulation chain coupled to the antenna and comprising 
multiple cascaded variable-delay polarization modulators, 
each of which may include a translation apparatus providing 60 
variable magnitude of phase delay between two orthogonal 
linear polarizations of an input signal, a beamsplitter and a 
first arm and a second arm optically coupled to one another 
via the beamsplitter. Each of the first and second arms may be 
formed of opposed reflective elements disposed on opposite 65 
sides of the beamsplitter, each of which may include at least 
one mirror coupled to the translation apparatus such that a 
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first distance separating the beamsplitter from at least one of 
the opposed reflective elements is variable. 

BRIEF DESCRIPTION OF THE DRAWING 

FIG. 1 is a simplified block diagram of an overview of a 
system configured to improve quasi -optical signal detection 
operations according to an embodiment of the invention. 

FIG. 2 is a simplified representation of a variable-delay 
polarization modulator element useful in the context of the 
system of FIG. 1 according to an embodiment. 

FIG. 3 is a simplified representation of a variable-delay 
polarization modulator element useful in the context of the 
system of FIG. 1 according to an embodiment. 

FIG. 4 is a simplified block diagram illustrating a variable- 
delay polarization modulator useful in the context of the 
system of FIG. 1 according to an embodiment. 

FIGS. 5 through 9 are graphical depictions of measured 
and simulated performance of a millimeter-wave or quasi- 
optical signal detection apparatus using the concepts devel- 
oped with respect to the variable-delay polarization modula- 
tor element(s) according to an embodiment. 

DETAILED DESCRIPTION 

In the following detailed description, reference is made to 
the accompanying drawings that form a part hereof, and in 
which are shown, by way of illustration, specific embodi- 
ments which may be practiced. These embodiments are 
described in sufficient detail to enable those skilled in the art 
to practice the embodiments, and it is to be understood that 
other embodiments may be utilized, and that logical, 
mechanical, electrical and other changes may be made, with- 
out departing from the scope of the embodiments. 

Ranges of parameter values described herein are under- 
stood to include all subranges falling therewithin. The follow- 
ing detailed description is, therefore, not to be taken in a 
limiting sense. 

The detailed description is divided into five sections. In the 
first section, a system level overview is provided. In the sec- 
ond section, a physical example of a variable-delay polariza- 
tion modulator element is presented. In the third section, tools 
for modeling variable-delay polarization modulator elements 
are developed. In the fourth section, hardware and an operat- 
ing environment in conjunction with which embodiments 
may be practiced are described. In the fifth section, a conclu- 
sion of the detailed description is provided. A technical effect 
of the systems and processes disclosed herein includes at least 
one of facilitating capability for measurement of polarization 
components of millimeter-wave or quasi-optical signals. 

§1. System Overview 

FIG. lisa simplified diagram of an overview of a modified 
system 100 configured to improve quasi-optical signal detec- 
tion operations. The system 100 may receive a signal 102 
from a source 103 via one or more antennae 105 (only one 
antenna is shown in FIG. 1), which may be parabolic, or may 
be any other type of antenna suitable for the wavelength of 
interest. A receiver 110 can be coupled to the antenna 105 and 
can provide received signals 102 via a signal port 115 to a first 
variable-delay polarization modulator element VPM 120, 
which may be coupled via signal path 122 to a second vari- 
able-delay polarization modulator element VPM 124, which, 
in turn, can pass the modulated signal on to other receiver 
elements (not shown in FIG. 1) via an output port 125. 
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It will be appreciated that while two variable-delay polar- 
ization modulator elements VPM 120, 124 are shown coupled 
in cascade in FIG. 1, additional variable-delay polarization 
modulator elements. In one embodiment, the system 100 may 
be included in the receiver 110. As well, one or more of the 
variable-delay polarization modulator elements may inte- 
grate into the system 100 in other ways, such as forming a 
portion of the antenna 105, or as coupling the antenna 105 to 
the receiver 110. Also, it will be appreciated that the receiver 
110 may be employed in capacities not associated with a 
signal reception system per se. For example, the cascaded 
variable-delay polarization modulator elements VPM 120, 
124 may be employed as a portion of a calibration apparatus 
or for other purposes. 

In one embodiment, the system 100 may be a radio recep- 
tion system configured to detect modulation present on a 
radio signal 102. Such modulation may result from a variety 
of sources, including man-made sources or naturally -occur- 
ring objects and associated phenomena. In other words, 
modulation present on the signals 102 may represent any of a 
variety of phenomena, including communications -related 
activities, remote sensing applications, naturally-occurring 
objects and associated phenomena and other scenarios 
wherein detection of such modulation is desirable. In one 
embodiment, the system 100 may be at least a portion of an 
astronomical instrument, such as a radio telescope, that may 
be employed for purposes of cosmology or other study. 

The variable-delay polarization modulator elements 120 
may stand in contrast to the conventional birefringent modu- 
lators in that the basis of phase introduction can be held fixed, 
but the magnitude of the delay may be variable. As a result, a 
variable-delay polarization modulator element 120 may pro- 
vide an adjustable delay between two orthogonal linear polar- 
izations of an electromagnetic wave comprising the signal 
102. There are many known elements in the art capable of 
providing such variable delay. For example, a Martin-Pup lett 
interferometer that does not include an input polarizer has 
been described in the past. 

In this invention, the polarization modulation can be 
explicitly separated from the polarized detection (obtained 
via an analyzer) of the signal 102. One advantage realized by 
this treatment may be that the basis of the variable-delay 
polarization modulator can be rotated at an arbitrary angle 
with respect to the analyzer. 

§11. Variable-Delay Polarization Modulator Elements 

FIG. 2 shows a simplified schematic diagram of a variable- 
delay polarization modulator 200 of a Martin- Pup lett inter- 
ferometer, i.e., an n=l case, for an angle of jt/4. The variable- 
delay polarization modulator 200 may include an input port 
202, and output port 204, a first rooftop mirror 205 compris- 
ing a first planar reflective surface 207 and a second planar 
reflective surface 209 at right angles to one another, a second 
rooftop mirror 210 comprising a first planar reflective surface 
212 and a second planar reflective surface 214 also at right 
angles to one another and a wire grid 215 acting as a beam- 
splitting reflector and as a polarizing filter interposed between 
the first 205 and second 210 rooftop mirrors. When d 1 =d 2 , the 
modulator 200 may behave as a mirror; when there is a path 
difference, it can change the polarization state of the incom- 
ing radiation 220. 

An incoming light beam 220 can enter through the input 
port 202 and then be incident on the grid 215. A portion 222 
of the beam may be transmitted through the grid 215 and may 
strike the first planar reflective surface 212 of the second 
rooftop mirror 210, giving rise to a first reflected beam 226 
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that may strike the second planar reflective surface 214 of the 
second rooftop mirror 210 and in turn gives rise to a second 
reflected beam 228 that may travel parallel to the transmitted 
beam portion 222, but in the opposite direction. The second 
5 reflected beam 228 may strike the grid 215, giving rise to a 
third reflected portion 230 that may travel parallel to the first 
reflected beam 226, but in the opposite direction, to the output 
port 204. 

Another portion of the incoming light beam 220 may be 
reflected by the grid 215, giving rise to a first reflected beam 
236 that can strike the first planar surface 207 of the first 
rooftop mirror 205, thus resulting in a second reflected beam 
238 that can strike the second planar surface of the first 
15 rooftop mirror 205, and in turn giving rise to a third reflected 
beam 240 traveling in the opposite direction. According to the 
example of an embodiment shown in FIG. 2, radiation may 
enter from the left (arrow 220), and may be split by a beam- 
splitting grid 215 into beams 222 and 236. These two com- 
20 ponents of polarization may then be directed to a respective 
one of rooftop mirrors 205 and 210, which can rotate the 
polarization by ninety degrees with respect to the wires form- 
ing the beam- splitting grid 215. 

Beams 230 and 240 can recombine at the beam-splitting 
grid 215 and may exit towards the top of the variable-delay 
polarization modulator 200. In the following analysis, the 
angle of the variable-delay polarization modulator 200 cor- 
responds to the angle of the beam-splitting grid 215. How- 
30 ever, one skilled in the art will recognize that other possibili- 
ties exist that may fall within the scope of this invention. 

The variable-delay polarization modulators 200 can be 
configured as follows: the variable-delay polarization modu- 
lator closest to the polarization- sensitive detector, i.e., vari- 
35 able-delay polarization modulator! , may have its beam-split- 
ting grid 215 oriented at forty -five degrees with respect to the 
axis of the detector (Q axis), and the second variable-delay 
polarization modulator, i.e., variable-delay polarization 
modulator 2 , may have its grid 215 oriented at an angle of 
40 twenty -two and a half degrees (i.e., jt/8) with respect to the 
detector axes. Full modulation of all linear and circular polar- 
izations states can be achieved with this configuration. 
Employing this architecture in a polarimeter that measures 
linear polarization may be modeled as described below. 
45 Denoting a first variable-delay polarization modulator, i.e., 
variable-delay polarization modulator! , which may be an 
input variable-delay polarization modulator and thus near the 
signal source, as having zero phase delay, and then switching 
a second variable-delay polarization modulator, i.e., variable- 
50 delay polarization modulator 2 , between delays of Oand jt, 
then the detector axes, projected onto the plane of the sky, can 
switch between Q and -Q. With variable-delay polarization 
modulator set to a phase delay of jt, switching variable-delay 
polarization modulator 2 between 0 and jt can project the 
55 detector axes to ±U. The dual variable-delay polarization 
modulators can provide two degrees of freedom, namely, the 
phase delay of each. The angles chosen for the two basis sets 
may be those for which the two degrees of freedom corre- 
spond to orthogonal coordinates on the Poincare sphere, thus 
60 allowing all polarization states to be accessible to the detec- 
tor. 

Qualities providing utility for this architecture as a candi- 
date technology for polarimeters in far-infrared through mil- 
limeter regions of the spectrum for astronomical observation 
65 applications may include: (i) this approach may allow imple- 
mentations that are able to the first reflected beam 236 and 
parallel to it. The third reflected beam 240 combines with the 
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third reflected beam portion 230 and the combined beams 23 0 
and 240 exit towards the top of the variable-delay polarization 
modulator 200. 

The polarizations of the electric E and magnetic H compo- 
nents of the various beam portions can be as noted in FIG. 2, 
and the rooftop mirrors 205 and 210 can act to rotate the 
polarization by ninety degrees with respect to the wires form- 
ing the beam-splitting grid 215. The path length of either the 
arm of the modulator 200 including the first rooftop mirror 
205 or the second rooftop mirror 210 may be modified, for 
example by translating the first rooftop mirror 205 parallel to 
the path noted by the arrow 250, resulting in modulation of 
relative phase of the third reflected beam 240 relative to the 
third reflected beam portion 230 by changing the path length 
d l5 or by translating the second rooftop mirror 210 parallel to 
the arrow 260, again resulting in controlled modulation of the 
relative phases of the third reflected beam portion 230 relative 
to the third reflected beam 24 0 by changing the path length d 2 . 
In either case, the third reflected beam portion 230 and the 
third reflected beam 240 may combine to provide either con- 
structive or destructive interference varying with the relative 
path length modulation of d 1 -d 2 . When d 1 =d 2 , the modulator 
200 behaves as a mirror; when there is a path difference, the 
polarization state of the incoming radiation 220 is modified as 
described infra. 

FIG. 3 is a simplified representation of a variable-delay 
polarization modulator element 300 useful in the context of 
the embodiment of FIG. 1. The variable-delay polarization 
modulator element 300 may include an input port 302, an 
output port 304, a mirror 305 and a polarizing wire grid 315 
that also acts as a partially reflective surface or beam splitter. 
An input signal 320 (solid line and arrow) may enter via the 
input port 302 and strike the grid 315, giving rise to a partially 
transmitted portion 322 (dashed line and arrow) and a par- 
tially reflected beam portion 330 (dot-dash line and arrow). 
The partially transmitted beam portion 322 may be incident 
on and reflected by the mirror 305, resulting in a reflected 
beam 340, at least a portion of which may traverse back 
through the grid 315 parallel to the reflected beam portion 330 
and which may combine with the reflected beam portion 330 
at the output port 304. Motion of the mirror 305 relative to the 
grid 315 or vice-versa may result in a variable path length 
difference for the beam portions 330 and 340, again giving 
rise to modulated constructive or destructive interference, as 
described above. 

In the following analysis, the angle of the variable-delay 
polarization modulator 200 may correspond to the angle of 
the beam-splitting grid 215. In contexts such as the receiver 
110 of FIG. 1, where multiple variable-delay polarization 
modulators 200 are coupled in cascade, the variable-delay 
polarization modulators 200 may be configured as follows: 
the variable-delay polarization modulator 200 closest to the 
polarization-sensitive detector, i.e., variable-delay polariza- 
tion modulator l3 may have its beam-splitting grid 215 ori- 
ented at forty-five degrees with respect to the axis of the 
detector (Q axis), and the second variable-delay polarization 
modulator 200, i.e., variable-delay polarization modulator 2 , 
may have its grid 215 oriented at an angle of twenty -two and 
a half degrees (i.e., Jt/8) with respect to the detector axes. Full 
modulation of all linear and circular polarizations states can 
be achieved with this configuration. Employing this architec- 
ture in a polarimeter that measures linear polarization may be 
modeled as described below. Denoting a first variable-delay 
polarization modulator, i.e., variable-delay polarization 
modulator!, which is an input variable-delay polarization 
modulator and thus near the signal source, as having zero 
phase delay, and then switching a second variable-delay 
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polarization modulator, i.e., variable-delay polarization 
modulator 2 , between delays of 0 and jt, then the detector axes, 
projected onto the plane of the sky, will switch between Q and 
-Q. With variable-delay polarization modulator! set to a 
5 phase delay of jt, switching variable-delay polarization 
modulator 2 between 0 and jt projects the detector axes to ±U. 
The dual variable-delay polarization modulators may provide 
two degrees of freedom, namely, the phase delay of each. The 
angles chosen for the two basis sets may be those for which 
to the two degrees of freedom correspond to orthogonal coordi- 
nates on the Poincare sphere, thus allowing all polarization 
states to be accessible to the detector. 

Qualities providing utility for this architecture as a candi- 
date technology for polarimeters in far-infrared through mil- 
15 limeter regions of the spectrum for astronomical observation 
applications include: (i) this approach allows implementa- 
tions that are able to cover the full Poincare sphere, in contrast 
to birefringent plates, which can assess either linear or circu- 
lar polarization, but not both; (ii) the variability of the path 
20 difference between orthogonal polarization states may facili- 
tate retuning for use at multiple wavelengths; (iii) frequency- 
dependent antireflective coatings can be avoided, in contrast 
to approaches relying on transmission through thick dielec- 
tric plates; (iv) the complexities of shafts, bearings and gears 
25 to rotate the birefringent plate in modulators based on such 
plates can be avoided. Each of these qualities may present a 
benefit in the context of measuring polarized flux of astro- 
nomical and cosmological sources from space-borne obser- 
vation platforms. 

30 The variable-delay polarization modulator in a Martin- 
Puplett interferometer is a specific example of more general 
cases. The Martin-Puplett interferometer is configured with 
the relative angle of the variable-delay polarization modula- 
tor at forty-five degrees with respect to the analyzer. The 
35 analytical expressions for the polarization components in this 
case may follow from those for the case of the single variable- 
delay polarization modulator placed at an arbitrary angle, 
although the physical implementation is different. 

The development and disclosure of aspects of expressions 
40 for the transfer function for multiple variable-delay polariza- 
tion modulator elements in cascade at arbitrary relative ori- 
entations may proceed from an example of two such ele- 
ments, as described in more detail below. Within that context, 
initially developing an expression for the transfer function for 
45 a single variable-delay polarization modulator then later may 
facilitate modeling of multiple cascaded variable-delay polar- 
ization modulators having arbitrary relative orientations. The 
discussion that follows assumes that the modulation passband 
is sufficiently narrow that the phase delays introduced at the 
50 center wavelength approximately apply over the whole band- 
width. 

§111. Modeling Variable-Delay Polarization 
Modulator Elements 

55 

A Mueller matrix representation (e.g., Eqs. (2) through 
(10) and Table I, infra) of the variable-delay polarization 
modulator 200, according to an embodiment, is now devel- 
oped in §III(a), based on consideration of the interior portion 
60 of a Martin-Puplett interferometer. Via this representation, 
which is given strictly by way of example and is therefore not 
limiting, the frequency-dependent performance of a variable- 
delay polarization modulator can be modeled in §III(b) (e.g., 
Eqs. (11) through (16), infra). An alternative embodiment of 
65 architecture for the variable-delay polarization modulator is 
described in §III(c), and systematics are considered in §III(d). 
Experimental results from laboratory tests of a single vari- 
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able-delay polarization modulator are then presented in §111 
(e) (and Eq. 17), and those empirical results are compared in 
§III(f) to simulations based on the transfer function models 
developed and disclosed herein (FIGS. 5 through 9). In §111 
(g), polarization matrix methods are described (e.g., Eqs. 18 5 
through 24 and Table II, infra). A summary of these aspects is 
provided in §III(h). 

A Martin-Puplett interferometer includes a variable-delay 
polarization modulator with an analyzer on the output end 
nominally oriented at an angle of forty -five degrees with 10 
respect to the beam-splitting grid 215. For spectrometer 
applications, one of the ports on the input side is shorted by a 
grid 215 oriented either parallel or perpendicular to the ana- 
lyzer (a grid 215 may have functionality other than beamsplit- 
ting). The model developed below of the interior of the Mar- 15 
tin-Puplett interferometer in terms of Jones and Mueller 
matrices allows a general variable-delay polarization modu- 
lator to be described at an arbitrary angle with respect to the 
optical system within which it is operative. The convenience 
of use of standard polarization matrices facilitates modeling 20 
of multiple variable-delay polarization modulators cascaded 
serially. 
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matrices for the radiation in each of the arms of the variable- 
delay polarization modulator 200, that is: 

J VPlJy^ 4) =j VPM ' 1 '’(jl/ 4) +J (2) 

In turn, each of the terms can be decomposed into a product 
of the Jones matrices of the elements in each optical path, for 
example, using the Jones matrices given below in Table I. A 
result of decomposing the terms of Eq. (2) as described is 
shown below in Eqs. (3) and (4): 

7vpm(*/ 4) = 7»T(>r/4)7My*r(0)7M)7i™(tf/4) = (3) 

(1 1 )exp(«4^i/A) 

Ul -lj 2 ’ 

and 

7$«w(*/4) = (4) 

_ _ _ _ _ f 1 -1 ^ exp(i4^i /A) 

J I ^)J z (d2)J z (d 2 )J 1 4) = i j ^ ■ 


TABLE I 


Summary of physical transformation models for optical elements, including Jones matrices. For 
linear distance transformation, d represents distance traveled; for the mirror, rotation has no 
effect; for the wire grid, 0 represents the angle of the grid wires with respect to the H-axis; 
for the rooftop mirror, 0 respresent the angle of the between the roofline and the H-axis; for the 
birefringent plate, 0 represents the angle between the fast axis of birefringence and the H-axis. 
and "E, is half of the phase delay introduced between the orthogonal polarizations. 


Description 

Symbol Jones matrix 

Stokes Expansion 

Linear distance 

J z (d) 

f exp(i2rrd/A) 0 'j 

, I exp(i2jtd/A) 



l o 

exp(i2rrd/A ) } 


Mirror 

Jm j 

f 1 °1 


Q 



lo lj 



Wire grid (ref.) 

JwR j 

( cos 2 

9 sin 0 cos 0 ) 

(Q + T cos20 + iV sin20)/2 



[ — sin 0 cos 0 sin 2 0 J 


Wire grid (trans.) 

J W T (0) j 

( sin 2 

0 — sin 0 cos 0^ 

(I + Q cos20 - U sin20)/2 



[ - sin 0 cos 9 cos 2 9 } 


Coord, rotation 

R(0) 

( cos 9 

sin0 ^ 

I cos0 +iV sin0 



l -sin0 

cos 9) 


Rooftop Mirror 

Jrt (Q) j 

( cos 2 9 

sin 20 \ 

I cos20 + iV sin20 



l -sin 29 

cos 20 J 


Biref. plate 

Jwp 1) j 

( cos - 

i sin £ cos 2 0 

-i sin f sin 2 0 A 



l -isin£sin20 -cos f -isin f cos 20 J 


Initially, a model is derived for the simple case of a vari- 
able-delay polarization modulator at an angle of forty-five 
degrees. This simple model is then generalized via similarity 60 
transformation. 

§III(a). Mueller Matrix Representation of a Single 
Stage 

65 

For the simple case, the Jones matrix J represent- 
ing this configuration can be expressed as a sum of Jones 


Defining As4jt(d 1 -d 2 )/X and setting 5= A/2, Eq. (5) fol- 
lows: 


J vpm (x / 4, f ) = 


exp(i27r(di +d 2 )/A)f cos(£) 
V isin(f ) 


-ficos(^) ' 

— COS (£) , ' 


(5) 


Eq. (5) forms a basis for derivation of an expression 
descriptive of a variable-delay polarization modulator 200 
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placed at an arbitrary angle 0. The definition of 0 describes the 
angle of the grid 215 with respect to H for the input radiation. 
To form this description, the coordinate system is trans- 
formed into the coordinate system used for the analysis 
above, the transformation for J FiW (jt/4) is applied, and then 
the coordinate system is transformed back. However, reflec- 
tions present a subtlety. For an odd number of reflections, the 
angle 0 of the device from the perspective of the outgoing 
beam is the negative of that from the perspective of the incom- 
ing beam. This aspect is included as shown below. Setting 
%=(0-jt/4), the description as stated in Eqs. (6) and (7) is 
provided: 


_ _ _ _ e\p(s27r(dt + d 2 ) / X) (6) 

Jvpuix* £) = R{-X)Jvpm{xI*)R{x) = 2 ~ 

( cos (if) + isin(f)cos(2^) -«sin(^)cos(2^) ) 

isin(^)cos(2^) -cos(£) + isin(£)sin(2*) / 

and 

Jvpm (0, £) = exp(i27r(di + d 2 ) / A) 0) 

( cos(£) - jjsin(f)cos(20) -isin(£)cos(20) A 
v isin(£)cos(20) -cos(£) - isin(£)sin(20) J 

In terms of the Stokes parameter basis set, this expression 
is as shown below in Eq. (8): 

7p7>A#(0» l)=exp (i 2it(d 1 +d 2 yk) \Q cos (l)-z' si n cos 

(20)+*Tsin(20))]. (8) 

Aside from a phase factor (irrelevant in measuring power), 
^vpm = Q^wp' Accordingly, the action of the variable-delay 
polarization modulator 200 is equivalent to that of a birefrin- 
gent plate (having its birefringent axis oriented at angle 0, 
with a delay A=2^) followed by a reflection (represented as 
the Jones matrix Q). 

The matrix in Eq. (8) is unitary and has a determinant of - 1 . 
Thus its Mueller representation is expected to describe the 
symmetries of the Poincare sphere. By expanding the density 
matrices in the Pauli matrix basis both before and after per- 
forming similarity transformations corresponding to the opti- 
cal system, the Mueller matrix for the system is determined, 
as shown below in Eq. (9): 


M vpm (Q, A) = 


(9) 


no o o \ 

0 cos 2 20 + cosAsin 2 20 -sin20cos20(l - cosA) sin20sinA 
0 sin20cos20(l - cos A) -sin 2 20 - cosAcos 2 20 -cos20sinA 
v 0 sin20sinA cos20sinA -cosA , 
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§III(b). Polarization Modulation 


It is assumed for purposes of this nonrestrictive example 
that the detector at the back end of the optical system is 
5 sensitive to Stokes Q. This is a statement about the orientation 
of the analyzer in the polarimetric system. Strictly speaking, 
a Q-sensitive detector requires a differencing of two orthogo- 
nal linearly polarized detectors with an orientation that we 
choose to define as the Q-axis. However, the model developed 
to below is applicable to the class of polarization detector strat- 
egies that collect only one linear polarization. Such detectors 
are technically sensitive to Q±I, but to lowest order, or in ideal 
modulation, I does not couple to the polarization modulation. 

The modulator changes the polarization state of the detec- 
15 tor as projected onto the sky. For a single variable-delay 
polarization modulator, the polarization state that the detector 
measured can be calculated from the second column of the 
Mueller matrix, as shown below in Eq. (11): 

Qdet = Qs/iy( cos2 2Q + cos A sin 2 20)+L r sA;v cos 20 sin 
AO 20(l-cos A)+r s ^(sin 20 sin A). (1 1) 

For a variable-delay polarization modulator, 0 is fixed and 
A is modulated. Using a single variable-delay polarization 
modulator, it is not possible to completely modulate Q, U and 
25 V. An example of this might be obtained by setting 0 to jt/4 . In 
this case, the expression of Eq. (12) is applicable: 


QderQsky cos A+ V sky sin A. (12) 

Here, Q and V are modulated, but U is not. This obtains as 
a result of the fact that at this grid 215 angle, U and -U 
propagate through the system separately, without interfering. 

An advantage of the variable-delay polarization modulator 
is that the phase freedom allows a straightforward method for 
calculating Q and V across large frequency bands. A conve- 
nient definition set is: A=k8, such that k=2jc/X, where 8 rep- 
resents the total path length difference between two orthogo- 
nal polarizations (for traditional Martin-PuplettPupplet beam 
paths 8=2(d 2 -d 1 )). The spectrum measured by the polarized 
detector is then dependent on 8 and is given by Eq. (13) 
below: 

Q\k, b)=Q(k)cos(kb)+V(k)s'm(kb). (13) 

For bolometric detectors, the signal is integrated over the 
passband of the instrument, 4>(k), as described in Eq. 14: 


45 


Q'(S)= Q'(k',6Mk')dk? 


(14) 


50 


[Q- (k')cos(k' 6) + V(k')sink'$]<p(k')dk' . 


The matrix shown in Eq. (9) may be expressed as a product of 
symmetry operations on the Poincare sphere, as given below 
inEq. (10): 

A) = YQvYQuRy('2B)RQ(b)Ry(-'2B) (10) 

= YqvYquM wp{0. A). 

Here, T QV represents a reflection about the Q-V plane, T QU 
T qv represents a reflection about the Q-U plane, R v repre- 
sents a rotation about the V-axis, and R e represents a rotation 
about the Q-axis. The matrix M vpm (0, A) is the Mueller 
matrix for a wave plate. 


Taking the Fourier transform of both sides of Eq. (14) 
yields Eq. (15): 

55 

r 6 2 ( 15 ) 

1 /2[Q(k) + iiV(k)]<f>(X) = 1/2 n\ Q'(d)e tk6 dd. 

60 The real part of the Fourier transform of the interferogram 
is the spectrum of Stokes Q from the source, while the imagi- 
nary part is the Stokes V spectrum. Note that broadband 
modulation relies on sampling a large enough range of path 
length differences. 

65 For implementation of a Fourier Transform Spectrometer 
via a Martin-Puplett interferometer, a horizontal or vertical 
grid 2 15 is placed at the input port of the device 20 0 . The input 
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polarization state is then purely Q, enabling the internal grid 
215 to function as a broadband, frequency-independent 
beamsplitter. This is equivalent to shorting one of the input 
ports. Thus, assuming that the polarization of the source is 
small. 


Q(k) = 1/2 (/(*)), 


and Eq. (15) reduces to the unpolarized spectrum of the 
source 103. In this case, the device 200 does not measure 
polarization, but relies on the fact that the input grid 215 is 
fixing the polarization to something that is known. 

A disadvantage to this architecture is its insensitivity to 
Stokes U. For experiments using space-borne measurement 
platforms, U can be recovered by rotation of the spacecraft. 
For experiments using terrestrial measurement platforms, 
sufficient rotation is problematic, and other techniques may 
be required. 

An alternative to measurement platform rotation is to cas- 
cade two instruments 200 serially. The functional form of the 
polarization signal may be found by simply chaining the two 
Mueller matrices together. 

The transfer equation of the system is now 
A 1 )M fpm (02, A 2 )S det , where modulator 2 is closer to the 
detector than modulator . Because the Q det detectors are sen- 
sitive only to Q, we solve for the second column of the 
resulting matrix, shown below in Eq. (16): 

&e? = 2s^[( cos2 20 i+cos sin 2 20 1 )(cos 2 20 2 +cos A 2 
sin 2 20 2 )-sin 20 x cos 20 x sin 20 2 cos 20 2 )(l-cos 
A 1 )(l-cos A 2 )+sin 20 j sin 20 2 sin Ai sin AJ+ 

U sJ -,[ sin 20 ! cos 20 1 (l _cos A 1 )(cos 2 20 2 +cos A 2 
sin 2 20 2 )-(sin 2 20 1 +cos A x cos 2 20 1 )sin 20 2 cos 
20 2 (1-cos A 2 )-cos 20 1 sin 20 2 sin A : sin A 2 J+ 

V sky [ sin 20 1 sin A 1 (cos 2 20 2 +cos A 2 sin 2 20 2 )+cos 
20j sin 20 2 cos 20 2 sin A^l-cos A 2 )-sin 20 2 cos 
AjSinA 2 ]. (16) 

A specific case, given by way of example, where B^Jt/8 
and 0 2 =0/4, is now considered. These angles are chosen to 
allow full sampling of the Poincare sphere. Polarized sensi- 
tivities for selected pairs of phase delay settings for the pair of 
modulators are shown in Table II. Here it is possible to fully 
characterize the polarization state, and a simple method for 
doing this is to adopt a single phase delay over the entire 
bandwidth. In this case, the modulators are set to the desired 
detector polarization sensitivity and a measurement is made. 
The measurement is repeated for each state, thus building up 
information about the polarization state of the source. 


TABLE II 


Mapping of Q det onto the sky, corresponding to selected values 
of Ai and At for dual modulators having; 0i = Jt/8 and 0o = ji/4. 

A: 

a 2 

Qdet 

0 

0 

Qsky 

0 

Jt 

“Qsky 

Jt 

0 

U sky 

Jt 

Jt 

-U sky 

0 

jt/2 

-V sky 

Jt/2 

0 

y2(Q^ y + u sky )+v sky A/2 

Jt/2 

Jt 

-^(Qdcy + UJ-V^ 

Jt 

jt/2 

V sky 

Jt/2 

jt/2 

(Qsky-U Sky )/J2 


One of the strengths of this modulator is an ability to 
modulate quickly between different polarization states, in 
turn providing a benefit by putting the polarization signal 
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above the 1/f knee of the instrument noise spectrum. It may 
also be possible to extend the bandwidth in a way similar to 
the single modulator above. As a result, it is possible to scan 
these modulators through a range of delays and thus extract 
5 the frequency-independent Stokes parameters. 

§III(c). Other Implementations 

The architecture exemplified via the illustration of FIG. 2 is 
10 not a unique implementation of a variable-delay polarization 
modulator. Several different arrangements of grids 215 and 
mirrors correspond to Jones matrices that differ only by an 
absolute phase from those that describe the Martin-Puplett 
interferometer as derived above. A simple example uses a 
15 polarizing grid 215 placed in front of a mirror. This design is 
similar in structure to a reflecting waveplate, but, in this 
design, modulation occurs via modulating the grid 2 15 -mir- 
ror spacing, rather than by spinning the device. This alterna- 
tive design for a polarizing interferometer has been previ- 
20 ously employed because it provides compact features and is 
relatively easy to construct. This implementation is useful in 
a dual modulator system application because it requires sig- 
nificantly less space in the optical path than many other sys- 
tems need. However, it does present several drawbacks, 
25 including at least: (i) an inability to achieve a zero path length 
condition results in a slight decrease in effective bandwidth; 
and (ii) the two polarizations on the output side are slightly 
displaced. This effect may be mitigated slightly by using the 
modulators at close to normal incidence and by placing the 
30 modulators as close as possible to an aperture/pupil. 

§III(d). Systematic s 

In developing a polarization modulator, one must consider 
35 the possibility of instrumental effects introduced by the 
action of the modulation. In a dielectric half-wave plate, such 
an effect arises from the absorption properties of a birefrin- 
gent material. Loss tangents for light polarized along the fast 
and slow axes are generally different. The result is a modu- 
40 lated signal that appears at twice the rotational frequency of 
the birefringent plate. For the dual variable-delay polarization 
modulator system, there are two important effects to consider, 
as follows: (i) different settings of the translational stage 
result in changes in edge illumination, potentially introducing 
45 a spurious polarization signal; and (ii) differential absorption 
by the grids 215 and the mirrors 207, 212 of the modulator 
200. One of these concerns may be obviated by restricting the 
use of such modulators to slow optical systems in which the 
beam growth through the modulator is minimal. Another of 
50 these concerns is mitigated via use of rooftop mirrors (e.g., 
rooftop mirrors 207, 212), where the angle of incidence is the 
same, for different modulator portions. Thus, the Fresnel 
coefficients for each of the two polarizations remains essen- 
tially constant during the modulation process. 

55 

§III(e). Experimental Results 

The description of experimental results comprises four 
parts. In part (i), the setup is described. In part (ii), experi- 
60 mental procedures are explained. In part (iii), experimental 
results are described. In part (iv), resonances are considered. 

(i) Setup 

65 FIG. 4 shows a simplified block diagram 400 illustrating a 
variable-delay polarization modulator 401 useful in the con- 
text of the system 100 of FIG. 1. The example of the embodi- 
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ment provided in FIG. 4 also is helpful in testing and evalu- 
ating performance of the polarization modulation of a test 
vehicle exemplifying the principles of operation of a Martin- 
Puplett interferometer in the context of the disclosed subject 
matter. 5 

In the block diagram 400 of FIG. 4, the modulator 401 is 
shown as having an input port 402 (having a reference plane 
labeled TJ, and an output port 403 (having a reference plane 
labeled T 2 ). The input port 402 is 302 may be coupled to a 
feedhom 404304, which, in turn, may provide a signal 405 to 
(analogous to the signal 102 of FIG. 1) to the modulator 402, 
where the signal 405 has polarization elements 406 (dashed 
line) and 407 (dot-dashed line) associated therewith. 

The beam 405 exiting the horn 404 attached to the input 
port one 402302 subsequently can be redirected via a reflec- 15 
tive element 412 and then can pass through a polarizing grid 
415 that has its wires oriented at a forty -five degree angle in 
projection. Each orthogonal polarization 406, 407 is 306, 307 
can then be launched down a separate arm of the device 
401301 and reflect off of one of two rooftop mirrors 420, 422, 20 
each of which rotates the polarization vector by ninety 
degrees, as shown above with reference to FIG. 2. The beams 
may recombine at the polarizer 415 and may be refocused via 
a second reflective element 425 into a second feed horn 
426326 which may be coupled to port two 404. 2 : " 

The experimental setup is shown as being symmetric, and 
so the reverse light path is identical. The rooftop mirrors 420, 
422 are 320, 322 may be placed on translational stages (not 
illustrated), such that their relative distance (i.e., d l and d 2 as 
noted in FIG. 2) may be adjusted. The frequency -dependent 
phase that corresponds to this path length difference is the 
parameter that determines the mapping between polarization 
states on either end of the device 401. Note that the system can 
support ~1 mode in each polarization. As a result, Gaussian 
analysis was used in designing the optics. 

From a microwave circuit perspective, this device 401 may 
be modeled as a four-port device, with the two ports 402, 404 
on either end of the device 401 being defined as the vertically 
and horizontally polarized electric field modes. In these 4Q 
experiments, a Hewlett-Packard HP8106D millimeter wave 
vector network analyzer (hereinafter “VNA”) (not illustrated) 
was used to measure scattering parameters between these 
modes. The calibration reference planes are shown (T x 1 and 
T 2 ) in FIG. 4. The VNA may be used to measure the two by 45 
two scattering matrices of pairs of these ports, so in order to 
reduce contamination of results, a respective orthomode 
transducer (hereinafter “OMT”) 432, 434 was placed at the 
back of each respective feedhorn 404, 426 and terminates the 
unused polarization with a matched load. 50 

For purpose of these measurements, it is useful to regard 
the end of the feedhorn 404 attached to one port of the VNA 
as the source (analogous to the source 103 of FIG. 1) and that 
attached to another port of the VNA as the detector. The 
polarization state of the source for the measured data (pre- 55 
sented infra with respect to FIGS. 5 through 9) was set to be 
vertically -polarized light (a pure Q state), by orienting the 
waveguide appropriately. On the detector side (i.e., corre- 
sponding to the second port 403), both V and H states were 
quantified in successive measurements by respectively omit- 60 
ting and adding a ninety-degree twist to the WR-10 
waveguide between the OMT 434 and the T 2 calibration point 
403. The loss of the twist 436 is measured as 0.15 dB. The 
calibrated difference between the power associated with H 
and V gives a measurement of Stokes Q at the detector. 65 
The bandwidth of the test setup is approximately 78-115 
GHz. At the low end of the band, the band edge is defined by 
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that of the W-band feed horns 404, 426, and at the high end, it 
is defined by return loss due to the OMTs 432, 434. 

The experimental data presented via FIGS. 5 through 9, as 
described in §III(f), infra, employ a pair of W-band feed horns 
404, 426 (25-27 dBi) that are collimated by ellipsoidal mir- 
rors 412, 425 (f=25 cm). Each of the two rooftop mirrors 420, 
422 reflect 320, 322 reflects a component of polarizations. 
The mirrors are mounted on transports employed to adjust the 
path lengths of the individual polarizations. The polarizing 
grid is mounted such that the wires make an angle of forty- 
five degrees with the roof lines of (i.e., lines between corre- 
sponding or opposed points on the surfaces of) the rooftop 
mirrors 412, 425 in projection. The dashed and dotted lines 
show positions of the beam radius (8.7 dB edge taper and 20 
dB edge taper, respectively) of a Gaussian beam propagating 
through the structure for a 26 dBi feed and a wavelength of 
three millimeters (1 00 GHz). A dotted outline illustrates loca- 
tion for a ninety degree twist 436 that converts the sensitivity 
of port two from V to H. 

(ii). Experimental Procedure 

The zero path length position is found by first measuring 
the signal in the V direction at an experimental condition 
where the S 21 parameter was flat across the band. A first null 
condition was then used to achieve fine adjustment. V and H 
were measured for twenty- seven combinations of positions of 
the two mirrors (i.e., such as the rooftop mirrors 207, 212 of 
FIG. 2 or 420, 422 of FIG. 4) having path differences (e.g., 
such as di-d-j) corresponding to twenty -four degree steps in 
phase for X=3 mm. FIGS. 5 through 8 show four measured 
spectra. These plots include the expected transmission spec- 
tra (H a(l-cos A) and V a(l+cos A)), adopting a global gain 
of 0.9 dB to account for the expected loss beyond the calibra- 
tion port. The return loss of the system is about twenty-six dB 
and can been seen in the H component in FIG. 4. The trans- 
mission efficiency of the horns 404, 426 used here is not 
constant across the band, and tends to roll off at low frequen- 
cies. 

(iii). Experimental Results 

The experimental setup is described mathematically by Eq. 
(12). In this case, 

Qdet = m A) - yV(A)/ (//(A) + fV( A)) (17) 

= QsourceCOsA + V so „ rce sillA, 


where A=4jt(d 2 -d 1 )/X. Here, H(A) and V(A) are the powers 
corresponding to S 21 when the twist is included and excluded, 
respectively. For each frequency, the relative gain factor, f, is 
calculated by fitting for the average values of the signals in the 
H and V configurations and then taking the ratio. 

For every frequency, it is possible to measure Stokes Q and 
V. The result of this analysis is shown in FIG. 7. We find that 
that average Stokes parameters measured over the 78-110 
GHz band are Q=-l .002±0.003 andV=0.001 ±0.013. There is 
some non-zero power in Stokes V near the high end of the 
band. It is unclear whether this is due to a systematic effect, or 
to an unknown source polarization. 



17 

(iv). Resonances 


US 7,412,175 B2 


In this setup, proper termination of the unused port at both 
the entrance (202, 402) and exit (204, 403) apertures is essen- 
tial, as even small reflections can introduce resonances. These 
resonances are an indication of level of uncertainty of phase 
control of the radiation propagating through the device 200, 
300, 400. In turn, such uncertainty leads directly to a fre- 
quency-dependent random mixing between the Q and V 
polarization states, and hence to a decrease in the precision of 
the Martin-Puplett interferometer as a polarimeter. 

To address such factors, it is noteworthy that systematic 
“noise” may be mitigated via various levels of termination of 
the unused polarization. For example, addition of the OMTs 
432, 434 in the signal chain reduced noise in the transmission 
from three dB to one dB. This setup also employed a conven- 
tional horizontal grid (not illustrated), analogous to grid 215, 
315, 415, at the mouth of the appropriate feedhom 404, 426 to 
redirect any residual H component to a conventional eccosorb 
beam dump (also not illustrated). Such reduced the noise in 
transmission to one- fourth dB and further reduced the aver- 
age coupling Q into V from four percent to under one percent. 

On a telescope, this problem is mollified by the fact that the 
source 103 component is nearly perfectly terminated in the 
sky. This greatly reduces phase uncertainties in the system, as 
well as need for excessive polarization filtering. 

An additional consideration is that the ability of these 
variable-delay polarization modulator devices 200, 300, 400 
to operate at room temperature also favors application as 
calibration devices. An input polarized signal is easily trans- 
formed to test polarization response of a precision polariza- 
tion sensor. These variable-delay polarization modulator 
devices 200, 300, 400 are capable of transforming an initially 
linearly polarized state into an elliptical polarization state. 

§III(f). Experimental Results vs. Simulations 

FIGS. 5 through 9 summarize experimental results 
obtained as described above with respect to § §(i) through (iv). 
The thick solid line in each of these graphical depictions 
represents the spectrum of the vertical (V) linear polarization 
as measured at port two of the vector network analyzer. The 
thin solid line accompanying this represents the spectrum of 
the corresponding horizontal (H) linear polarization, also 
measured at port two of the vector network analyzer. The H 
polarization component is measured by adding a ninety -de- 
gree twist in the WR-1 0 waveguide attached to port two of the 
VNA. Predictions for H and V provided via the theory and 
models discussed above are plotted as thick and thin dashed 
lines, respectively. FIGS. 5 through 8 employ axes calibrated 
in decibels (dB, ordinate), and frequency, expressed in giga- 
Hertz (abscissa). 

FIG. 5 provides a graph 500 of signal transmission ampli- 
tudes in comparison to respective predicted expectation val- 
ues 502, 504 for horizontal (H) 512 and vertical (V) 514 
polarization components of a signal (e.g., signal 102 of FIG. 
1). The relative separation d : -d 2 of the rooftop mirrors 420, 
422 (FIG. 4) for the data displayed in FIG. 5 is a negative 
thirteen micrometers. 

FIG. 6 provides a graph 600 of signal transmission ampli- 
tudes in comparison to estimated values 602, 604 for respec- 
tive horizontal (H) 612 and vertical (V) 614 polarization 
components of the signal of FIG. 6. The relative separation 
d ] -d 2 of the rooftop mirrors 420, 422 (FIG. 4) for the data 
displayed in FIG. 6 is five hundred thirteen micrometers. 

FIG. 7 provides a graph 700 of signal transmission ampli- 
tudes in comparison to estimated values 702, 704 for respec- 
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five horizontal (H) 712 and vertical (V) 714 polarization 
components of the signal of FIG. 7. The relative separation 
d i -d 2 of the rooftop mirrors 420, 422 (FIG. 4) for the data 
displayed in FIG. 7 is a negative one thousand thirteen 
5 micrometers. 

FIG. 8 provides a graph 800 of normalized Stokes param- 
eters Q 802, 804 for respective horizontal (H) 812 and vertical 
(V) 814 polarization components of the signal of FIG. 8. The 
relative separation d 1 -d 2 of the rooftop mirrors 420, 422 
to (FIG. 4) for the data displayed in FIG. 8 is a negative one 
thousand thirteen micrometers. 

FIG. 9 provides a graph 900 of normalized Stokes param- 
eters Q 902 and V 904 as a function of frequency, by fitting to 
the twenty-seven positions for the rooftop mirrors 420, 422 of 
15 FIG. 4. Mean values for Q 902 and V 904 across the seventy - 
eight to one hundred fifteen GHz band are -1 .002 ±0.003 and 
0.001 ±0.01 3, respectively. 

§III(g). Polarization Matrix Methods 

20 

The overview of polarization matrix methods is developed 
in three parts. In part (i), Jones matrices are described. In part 
(ii), density matrices are described. In part (iii), Mueller 
matrices are described. 

25 

(i). Jones Matrices 

Jones matrices provide convenient modeling tools for 
analysis of radiation as it propagates through an optical sys- 
30 tern in which knowledge of phase is relevant. For an ideal 
case, it is assumed that all ports are impedance matched, and 
so no cavities are formed. This fonnulation is applicable for 
coherent radiation; however, it may be extended using the 
closely-related formalism of density matrices to treat situa- 
35 tions involving partially-polarized radiation or light. 

Jones matrices are two by two matrices expressing infor- 
mation regarding how orthogonal electrical field components 
transform in an optical system. The input Jones vector is 
expressed as shown below in Eq. (18): 

40 
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The output vector from an optical system may then be repre- 
sented by IEy>=JIEj>, where J is the vector transformation 
introduced by the optical system. The power measured at a 
detector at the back end of such a system may be expressed as 
50 shown in Eq. (19), below: 

<£ y j7^l£ / >=<£ 2 -l7 f '7^£ z >. (19) 

The matrix l det depends on properties of the detector used to 
55 make the measurement. In Jones matrix representation, 
Stokes parameters are represented by appropriate Pauli matri- 
ces and the identity matrix, as shown below: 
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-continued 


In these representations, the bar over the Stokes symbol indi- 
cates the Jones matrix representation. An un-barred Stokes 
symbol represents measurable power (e.g., Q=<EylJ^IEy>). 
The measured power in each of the Stokes parameters is given 1 0 
via Eqs. (21) through (24) below: 


U = (E\U\E) = 2 s Ji (E* h E v ); 
V = (E\V\E) = 23 (E* h E v ). 


Eqs. (21) through (24) connect the Jones matrix formula- 
tion of the Stokes parameters to their familiar definitions. 
These four Stokes matrices have the following multiplicative 
properties: if ( a q, a[, a^, UVV), a 0 a a =o^o^=a^ 25 

for a€(0, 1 , 2, 3) and and for j, k, 1 e(l, 2, 

3). These four matrices form a convenient basis for expressing 
Jones matrices. 

. . 30 

(n). Density Matrices 

Density matrices are two by two matrices which fully 
characterize the polarization state of the light. For the general 
case of partially polarized light, polarization arises because of 3 5 
time-averaged (statistical) correlations between the electric 
field components. A density matrix D is represented as shown 
below with reference to Eq. (25): 


(E* X E X ) (E* x E y ) 
(E;E Z ) (E* y E y ) 


The brackets < > indicate a time average. When the density 45 
matrix is expressed as a linear combination of Pauli matrices, 
D=Ia 0 +Qa 1 +Ua 2 +Va 3 , the coefficients are the Stokes 
parameters. 

Transformation of the polarization state by an optical sys- 
tem may be described via a similarity transformation, 1 )'=T^D 50 
J. Here, J is the Jones matrix describing the optical system. 

For the purposes of this disclosure, the manner in which the 
polarization state of the detectors maps onto the sky is rel- 
evant. Accordingly, D sky =VD det J. The transformation of the 
density matrix D and the expression for total power in the 55 
Jones formalism as expressed in Eq. (19) are notably similar. 

(iii). Mueller Matrices 

Mueller matrices are four by four matrices . In an analogy to 60 
special relativity, the inhomogeneous Lorentz group may be 
represented by a group of four by four real matrices acting on 
a Stoke vector, S=(I, Q, U, V), and these matrices are known 
as Mueller matrices. 

In the prior discussion, no limitations had been placed on J, 65 
the matrix describing the optical system being modeled. 
When the magnitude of the determinant of J is unity, there is 


a homomorphism between the group of two by two matrices 
having IdetJI T and the Poincare or inhomogeneous Lorentz 
group. With respect to the subject matter of this disclosure, 
the Mueller matrix that maps the Stokes parameters at the 
detector to the sky are of note, viz. , S sky =MS det . 

For polarization modulation, the instances for which the 
Jones matrices describing the optical system are unitary are of 
interest. For this case, Stokes I decouples from the other 
Stokes parameters, and the quantity P 2 =Q 2 +U 2 +V 2 is pre- 
served. This subgroup may be represented by three by three 
orthogonal submatrices representing symmetries on the sur- 
face of a sphere in space having Stoke I, Q and V as axes. This 
sphere is the Poincare sphere. 

When the group of density matrices are restricted to those 
with positive determinants, the system is described by SU(2), 
and so there is homomorphism between this group and the 
SO(3) group of rotations on the Poincare sphere. These are the 
groups that are important to a wave plate; however, the physi- 
cal reflection involved in the VPM architecture introduces a 
negative determinant, resulting in combinations of rotations 
on the Poincare sphere. 

§III(h). Summary 

Techniques for polarization modulation in which n phase 
delays between linear orthogonal polarizations are placed in 
series with arbitrary relative orientations as are described and 
modeled in this disclosure. The n=l and n=2 cases are spe- 
cifically considered, and it is noted that for appropriate rela- 
tive orientations, it is possible to fully modulate the polariza- 
tion in the n=2 case. In the far-infrared through millimeter 
where bandpasses are typically AA7L~0. 1 , this device can be 
used in a similar manner to a half-wave plate. Broader pass- 
bands (AX/X-0.3) may be accommodated using more com- 
plex modulation schemes. This architecture enables construc- 
tion of a modulator that may be made robust, broadband, and 
easily tunable to different wavelengths to provide frequency 
diversity. In addition, it permits complete determination of 
the polarization state of the incoming radiation by measure- 
ment of Stokes Q, U and V. 

§V. Conclusion 

A polarization sensitive receiver system is described. 
Although specific embodiments have been illustrated and 
described herein, it will be appreciated by those of ordinary 
skill in the art that any arrangement which is calculated to 
achieve the same purpose may be substituted for the specific 
embodiments shown. This disclosure is intended to cover any 
adaptations or variations. For example, although described in 
procedural terms, one of ordinary skill in the art will appre- 
ciate that implementations can be made in a procedural 
design environment or any other design environment that 
provides the required relationships. 

In particular, one of skill in the art will readily appreciate 
that the names or labels of the elements and apparatus are not 
intended to limit embodiments. Furthermore, additional pro- 
cesses and/or apparatus can be added to the components, 
functions can be rearranged among the components, and new 
components to correspond to future enhancements and physi- 
cal devices used in embodiments can be introduced without 
departing from the scope of embodiments. 
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What is claimed is: 

1. A millimeter wave receiver comprising: 

an input antenna; 

a polarization modulation chain coupled to the input 5 
antenna and comprising a plurality of cascaded variable- 
delay polarization modulators each including 

an input having a fixed basis of phase introduction; 

a translation apparatus providing variable magnitude of 10 
phase delay between two orthogonal linear polarizations 
of an input signal; 

a beamsplitter; and 

a first arm and a second arm optically coupled to one 
another via the beamsplitter, each of the first and second 
arms being formed of opposed reflective elements dis- 
posed on opposite sides of the beamsplitter, the opposed 
reflective elements in each arm including at least one 
mirror coupled to the translation apparatus such that a 20 
first distance separating the beamsplitter from at least 
one of the opposed reflective elements is variable. 

2. The receiver of claim 1, wherein one of each of the 
opposed reflective elements provides a collimation function. 25 

3. The receiver of claim 1, wherein one of each of the 
opposed reflective elements is formed as an ellipsoidal mirror 
and provides a collimation function. 

4. The receiver of claim 1, wherein one of each of the 
opposed reflective elements includes several planar surfaces 
oriented at right angles in a rooftop configuration and realiz- 
ing a reflected light path parallel to a path of incident light. 

5. The receiver of claim 1, wherein the beamsplitter 

includes a wire grid. 35 

6. The receiver of claim 1, wherein one of each of the 
opposed reflective elements is formed as an ellipsoidal mirror 
and provides a collimation function, another of each of the 
opposed reflective elements includes several planar surfaces 40 
oriented at right angles in a rooftop configuration and realiz- 
ing a reflected light path parallel to a path of incident light, 
and wherein the beamsplitter includes a wire grid. 

7. The receiver of claim 1, wherein the beamsplitter also 

functions as a polarizing filter. 45 

8 . The receiver of claim 1, wherein the variable-delay 

polarization modulators are capable of assessing both circular 
and linear polarization and provide capability for use at mul- 
tiple wavelengths. 50 

9. A method of characterizing the polarization of a signal 
comprising: 
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providing an input antenna; 

coupling a polarization modulation chain to the input 
antenna; said polarization modulation chain comprising 
a plurality of cascaded variable-delay polarization 
modulators each including: 
an input having a fixed basis of phase introduction; 
a translation apparatus providing variable magnitude of 
phase delay between two orthogonal linear polarizations 
of an input signal; 
a beamsplitter; and 

a first arm and a second arm optically coupled to one 
another via the beamsplitter, each of the first and second 
arms being formed of opposed reflective elements dis- 
posed on opposite sides of the beamsplitter, the opposed 
reflective elements in each arm including at least one 
mirror coupled to the translation apparatus such that a 
first distance separating the beamsplitter from at least 
one of the opposed reflective elements is variable; 
receiving a signal via the input antenna; and 
modulating the signal through the polarization modulation 
chain to characterize the polarization of the signal. 

10. The method of claim 9, wherein one of each of the 
opposed reflective elements provides a collimation function. 

11. The method of claim 9, wherein one of each of the 
opposed reflective elements is formed as an ellipsoidal mirror 
and provides a collimation function. 

12. The method of claim 9, wherein one of each of the 
opposed reflective elements includes several planar surfaces 
oriented at right angles in a rooftop configuration and realiz- 
ing a reflected light path parallel to a path of incident light. 

13. The method of claim 9, wherein the beamsplitter 
includes a wire grid. 

14. The method of claim 9, wherein one of each of the 
opposed reflective elements is formed as an ellipsoidal mirror 
and provides a collimation function, another of each of the 
opposed reflective elements includes several planar surfaces 
oriented at right angles in a rooftop configuration and realiz- 
ing a reflected light path parallel to a path of incident light, 
and wherein the beamsplitter includes a wire grid. 

15. The method of claim 9, wherein the beamsplitter also 
functions as a polarizing filter. 

16. The method of claim 9, wherein the variable-delay 
polarization modulators are capable of assessing both circular 
and linear polarization and provide capability for use at mul- 
tiple wavelengths. 



